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Abstract: We complete the Hamiltonian analysis of specific model of non- linear massive 
gravity that was started in arXiv:1112.5267. We identify the primary constraint and cor- 
responding secondary constraint. We show that they are the second class constraints and 
hence they lead to the elimination of the additional scalar mode. We also find that the 
remaining constraints are the first class constraints with the structure that corresponds 
to the manifestly diffeomorphism invariant theory. Finally we determine the number of 
physical degrees of freedom and we show that it corresponds to the number of physical 
modes of massive gravity. 
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1. Introduction and Summary 

Recently de Rham and Gobadadze proposed in Q an interesting formulation of the massive 
gravity wtiich is ghost free in the decoupling limit. Then it was shown in |Q] that this action 
that was written in the perturbative form can be resumed into fully non-linear action ^. 
The general analysis of the constraints of given theory has been performed in [^. It was 
argued there that it is possible to perform such a redefinition of the shift function so that 
the resulting theory still contains the Hamiltonian constraint. Then it was argued that 
the presence of this constraint allows to eliminate the scalar mode and hence the resulting 



theory is the ghost free massive gravity. However this analysis was questioned in |31| where 
it was argued that it is possible that this constraint is the second class constraint so that the 
phase space of given theory would be odd dimensional. On the other hand in the paper |32| 



very nice analysis of the Hamiltonian formulation of the most general gauge fixed non-linear 
massive gravity actions was performed with an important conclusion that the Hamiltonian 
constraints has zero Poisson brackets. Then the requirement of the preservation of this 
constraint during the time evolution of the system implies an additional constraint. As a 
result given theory has the right number of constraints for the construction of non-linear 
massive gravity without additional scalar mode ^. 

The Hamiltonian analysis of the manifestly diffeomorphism invariant non-linear mas- 



sive gravity with Stiickelberg fields was performed in [29| where corresponding Hamiltonian 
was found. Then using the observation firstly published in it was shown that this the- 
ory possesses one primary constraint. Unfortunately the presence of this constraint makes 
the calculation of the Poisson brackets between constraints very difficult due to the ab- 
sence of the inverse of the matrix V^^ = g^Wi4>^dj(j)^ and we were not able to perform 
this analysis for the case of four dimensional non-linear massive gravity. On the other hand 
such an analysis was performed for the case of two dimensional non-linear massive gravity 



iPor related works, see |, g g 1 1, |l§, 0, 0, g[ 0, [T8[ 0, |2l| g ||, p5| , ||, ||, @. 

^Alternative arguments for the existence of an additional constraints were given in |p3| even if the 
Hamiltonian analysis was not complete and the minimal non-linear massive gravity action was considered 
only. 
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with conclusions that there are no physical degrees of freedom left with agreement with 

II, mi- 



In [41 1 we analyzed the model of the non-linear massive gravity action introduced in 
1 33] written in the Stiickelberg formalism. This analysis was then reconsidered in with 
conclusion that this theory is free from the ghosts. 



The goal of this paper is to complete the analysis of the non-linear massive gravity 
action presented in [^]. We find the Hamiltonian for given theory and identify primary 
constraints. Then we rewrite the Hamiltonian to such a form where the scalar part of the 
Hamiltonian constraint will be proportional to the trace of the square root of the regular 
matrix. Then it would be possible to use the standard formula for the variation of the 
trace of the square root of the regular matrix and calculate corresponding Poisson brack- 
ets. Then we can analyze the requirement of the preservation of the primary constraints 
during the time evolution of the system and hence identify corresponding secondary con- 
straints. Finally we will check the stability of all constraints during the time evolution 
of the system. We find that the Hamiltonian and diff'eomorphism constraints are still the 
first class constraints and obey the basic principles of geometro dynamics [36, 38, |39|| . 
On the other hand we show that the additional primary constraint together with corre- 
sponding secondary constraint are the second class constraints and that these constrains 
could eliminate one additional degree of freedom so that the number of physical degrees of 
freedom correspond to the case of the massive gravity. In other words our results are in full 
agreement with the conclusion presented in j^^. However we mean that result derived in 
this paper is non-trivial and should be considered as an independent check of the absence 
of the ghosts in given theory due to the fact that we analyze theory without additional 
auxiliary fields so that the Hamiltonian analysis presented her is different from the analysis 
presented in 



We should also stress that our treatment has its own limitation due to the fact that 
we restrict to the case of one specific model of non-linear massive gravity action. It turns 
out that the extension of given analysis to the more general form of the non-linear massive 
gravity actions is very difficult due to the complicated relation between canonical momenta 
and time derivatives of the scalar fields. Unfortunately we were not able to find an inverse 
mapping that would allow us to write the Hamiltonian as a function canonical variables 
in these cases. It would be very interesting to find such Hamiltonian formulation and 
corresponding primary constraints between Stiickelberg fields even for the most general 
form of the non-linear massive gravity action. We hope to return to this problem in future. 



The structure of given note is as follows. In the next section (Q) we review some basic 
facts about the non- linear massive gravity action in the formulation presented in [ 2£ ] . Then 
in section (^) which is the main body of this paper we perform corresponding Hamiltonian 
analysis. We also identify primary and the secondary constraints and determine the number 
of the physical degrees of freedom. 
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2. Non-Linear Massive Gravity 

Our goal is to study non-linear massive gravity action in the following form ^ 

S = Ml j d?^dtN^ KijO'^^^Kki + ^^^R-m^TiA^ , (2.1) 
where we used 3 + 1 notation |34] ^ and write the four dimensional metric components as 



goo = -N' + Nig'^Nj , got = Ni , gij = gij , 

(2.2) 

Note that ^^^R is three-dimensional spatial curvature, Kij is extrinsic curvature defined as 

Kij = ^(dtgij - ViNj - VjNi) , (2.3) 

where Vj is covariant derivative built from the metric components gij and Q^^'^^ is de Witt 
metric defined as 

g^Jkl^hg^kgjl^g^lgJk^_g^Jgkl ^^^^^^^ 



^ ' ' - -iklmn ^ 



with inverse 

Qijki = -^(dikgji + gugjk) - -^gijgki , GijkiG"''"''' = -^(^T'Sj' + Si 5]^) . (2.5) 
Finally note that the matrix is defined as 

A^B = -V„0^V„(^B + g'^di(t>%(t>B , V„(/<^ = ^{M^ - N'dict)^) (2.6) 



and the trace defined in (2J.) is the trace over Lorentz indices A,B,C, . . . = 0, 1, 2, 3. 

We see that the action contains the potential term which is the square root of the 
matrix which can be defined as 

{VA)^s{^fc = ^^c ■ (2.7) 

For further purposes it is crucial to presume that A is regular matrix. Then when we 
perform the variation of this expression and multiply by (\/A)~^ from the right we obtain 

5{VA)^^ + iVA)^cS{^fDa^)-TB = S{A)^c{i^rTB ■ (2.8) 
Taking the trace the equation (|2.8| ) we immediately obtain ^ 



<5TYiVA = i(5A^s((VA)-i)V (2.12) 

This is the key formula that is used in the calculation of the Poisson brackets as we will 
see in the next section. 



^We use notations introduced in the paper |2£| ] 
''For review, see ISq]. 



^Note also that due to the matrix nature of objects A and B the following relation is not valid 

\/AB = VaVB (2.9) 
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3. Hamiltonian Analysis 

In this section we perform the Hamiltonian analysis of the action ( |2.1|) . For the General 
Relativity part of the action the procedure is standard. Explicitly, the momenta conjugate 
to A^, are the primary constraints of the theory 

7riv(x) « , 7ri(x) « (3.1) 

while the Hamiltonian takes the form 

(3.2) 

where vr*-? are momenta conjugate to gij with following non-zero Poisson brackets 

{5i,(x), 7r'='(y)} = i (5f 4 + 6i6^) 5(x - y) . (3.3) 

Note also that ttn^t^i have following Poisson brackets with N^N"^ 

{iV(x), ^^(y)} = 5(x - y) , {7V^(x), 7r,(y)} = ^^^(x - y) . (3.4) 

Now we proceed to the Hamiltonian analysis of the scalar field part of the action. Note 
that in 3 + 1 formalism the matrix A"^ takes the form 

A^B = -Vn<t>^yn't>B + g'^di^'^d.^B = K\ + vj,, (3.5) 

where 

K\ = -V„</)^V„0B , Kab = VAC'K^B = Kba , 



BA 



(3.6) 



Then the conjugate momenta pA are equal to 



PA - V^Tl^lA j (J 



^V^(Vn<^c(A-V2)C + ^^^(A-V2)C v.cA^) , = (Va)-i . 

(3.7) 



unless A and B commute. On the other hand since obviously A and A ^ commute the equation ( |2.9| ) gives 

VaVA-1 = I (2.10) 

which implies following important relation 



Va) = VA-i . (2.11) 
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Note that using the symmetry of Aab = A ba we can write (^^) in simpler form 

PA = M2m2V^(A^i/2)^BV„0^ . (3.8) 
Using this expression we derive following relation 

= {A-y^UciV^^ - A^^){A~y')nB 

which implies 



(3.9) 



Uab = {A-'/')AcV''''{A-'/^)nB , (3.10) 
where we introduced the matrix Uab defined as 

^AB = IjJlPAPB + VAB ■ (3.11) 



Note that when we multiply ( 3.10| ) by V from the right we obtain (we use matrix notation) 



UV = (A-i/V)(A-i/V) (3.12) 

which implies 

A-i/V = VlW . (3.13) 

This relation will be important below. The crucial point for the Hamiltonian analysis of 
the non-linear massive gravity is the fact that V^^ has the rank 3 as was firstly explicitly 
stressed in |^^. In fact, if we introduce the 4x3 matrix = di<f>^ and its transpose 
matrix (W'^Y = dicf)^ which is 3 x 4 matrix we can write 

V^^ = Wjg'^{W^)j'^ . (3.14) 

Then since W'l, g^^ have the rank 3 we obtain that V^^ has the rank 3 as well. As a result 
det V = 0. In other words V is not invertible matrix. 

With the help of these results it is easy to determine corresponding Hamiltonian 

= M'^PA - Csc = M^m^^NV^''{A~'/^)BA + N'pAd,(^^ = 



(3.15) 



usmg ( P^ and using an obvious relation Tr^-v/FII = Tvl^IW . With the help of these 



results we find the Hamiltonian for the action ( ^.1[) in the form 

H = j d^x{NnT + N'Tii + v'lTi + V^TTN + VcC) , (3.16) 

where 

-Ht = n^^ + , -Hi = nf^ + nf , (3.17) 
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and where tTj ~ , tttv ~ are the primary constraints of the theory. Note also that the 
Hamiltonian (3.16) contains primary constraint C whose exphcit form fohows from ( 3.1C| ) 
when we calculate the determinant of the matrix Hab- Using 



detn^B = -(l + ^^pV 



(3.18) 



and using ( |3.1[1| ) together with the fact that det V = we derive primary constraint C in 
the form 



C : 1 + -V^PAP^ « . 



(3.19) 



It is also important to stress that using the definition of Hab and the existence of the 
constraint C we obtain an important relation 



P^^AB = ( nrraPV + l]pB=CpB^O. (3.20) 



1 

Now we analyze the requirement of the preservation of the primary constraints. As usual 
the requirement of the preservation of the primary constraints VTiv ~ ,7ri w implies an 
existence of the secondary constraints 

nrix) « , Tiii^) ^ . (3.21) 

For further purposes we introduce the smeared form of these constraints ( p. 211 ) 

Tt{N) = J d^xNTir , Ts{N') = j d^xN'Tii . (3.22) 

It is not easy to determine the time evolution of the constraint C due to the fact that HV 
is singular matrix. To proceed let us express the trace of the square root of the matrix as 
power series in the form 

TtlVuV = ^ CnTvLil^Vr . (3.23) 
n 

Now we can write 

TtlUV = Yl^''di(t>''g'^d,(t>A = djct^AH^^'dicjiBg'' = H/ = Tr,n , 



where the trace Tr^ is the trace over spatial indices i,j,k... 
of (|3.24| ) it is easy to see that 



(3.24) 

1,2,3. Then with the help 
(3.25) 
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Now n is 3 X 3 matrix with the rank equal to 3 which imphes an existence of the inverse 
matrix 11"^. As a result we can easily determine the variation of the trace of the square 
root of given matrix 

STisVfl = ^TrsSflVfl-^ . (3.26) 
Then we can determine following Poisson brackets 



/ '^TrVn(y) l^lV(y) /— 



(x) 2 5(/)^(x) 



(3.27) 



= 2 m4M4 i^^^^PKdkcl^'' + ^»'/'^Pi^^^fc<^^)/^Vn-Hy))/'^(x - y) . 
Using these results we find 

{T5(iV^),C(x)} = -iV'5,C(x) 
and also 

{TT(Af),C(x)} = --g^p-' (8i[]VV9wn™8tfe9''-'V5^/l - 

(3.29) 



(3.28) 



where we used ( |3.20| ) and where we defined C^^ as 

, i« 2 



(3.30) 



where we defined ^vfr^j = yyfi^j = vft^ fcfi''^"'- Now it is easy to see that the 
requirement of the preservation of the constraint C during the time evolution of the system 
implies following secondary constraint 

dtC = {H,C} « {Tt{N),C} = v#^C^' « . (3.31) 

In summary, the theory possesses following collection of the primary constraints ttat w 
, TTj , C and secondary constraints Ht ~ « ,C^^ ^ 0. As a result the 
total Hamiltonian has the form 

Ht = J d^^iNUT + N'Ui + VNT^N + v'T^i + vcC + VHi + TcC^^) , (3.32) 
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where vn,v^,vc,T^,Tc are corresponding Lagrange multipliers. 

As the final step we have to analyze the preservation of all constraints. Note that in 
case of the General Relativity part of the constraints we have following Poisson brackets 



^GK(x)^<^(x-y)-^bH(y)|T5(x-y) 
d 



(3.33) 



The calculation of the Poisson brackets that contains scalar phase space degrees of freedom 
is more involved. However it is easy to find the Poisson bracket between generators of 
spatial diffeomorphisms 



{?^r(x),^f(y)} 



(3.34) 



that together with the Poisson bracket on the third line in ( 3.331 ) implies following form of 
Poisson bracket between smeared form of the diffeomorphism constraints 



It is also easy to see that 

{Ts{N'),Tt{N)} =TT{dkNN'') . 



(3.35) 



(3.36) 



Now we proceed to the calculation of the Poisson bracket {Tr(A^), Tr(M)}. By definition 
we have 



{Tf(iV),Tff(M)}= / d3xdViV(x){?^f(x),-?^ff(y)}M(y) 



M^m* / d^:>^d^yd'^zN{:>^)M{y) I ^(x) 



5(Vn)/(x)'5(vn)/(y) 



6px (z) 



(z) 



-Vaiy) 



(5(Vn)/(y)5(Vn).-'(x) 



5px{z) 6cl)^iz) 



where we used ( 3.27 ). This result together with (|3.33| ) implies ^ 

{Tt{N), Tt{M)} = TsiiNdjM - Mdj)g^') 



(3.37) 



(3.38) 



®It is important to stress that {Tg^(iV), T5?(M)} + {T??(A/), Tg^(iV)} = due to the fact that T?? 
does not depend on the spatial derivatives of Qij . 
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Now we have to determine whether ah constraints are preserved during the time evolution 
of the system. Let us now start with the primary constraints 7r^,Tri,C ~ 0. The case of 
the constraints vttv ~ , tTj ~ is trivial. Further, the requirement of the preservation of 
the constraint C gives 



dtCi^) = {i^T,C(x)} « J d\ ({iV(y)?^T(y),C(x)} + rc(y) {C^^(y), C(x)}) 

= M>^^"W + /^Vrc(y){c-(y),c(x)}. 



I d3yrc(y){c^'(y),c(x)}=o 

(3.39) 



that has the solution Tq = using the fact that 



{C^^(x),C(y)} = 2p^(x)a.n^^(x)V^(x) [VfF^'' (x)9,-5(x - y)j^^^^ +■■■ 

(3.40) 

does not vanish on the constraint surface. Note that . . . mean additional terms that arise 
from the explicit calculations of given Poisson brackets. In other words C and C^^ are the 
second class constraints. 

Now we come to the requirement of the preservation of the secondary constraints. Let 
us begin with the diffeomorphism constrains Tii or their smeared forms. Since C^^ ~ 
is manifestly diffeomorphism invariant we have {Ts(A^*), C^^} ~ and also using ( 3.2^ ) 
together with ( 3.35| ) and ( p. 36 ) we find that Hi is preserved during the time evolution of 
the system. In case of T-Lt we find that its time development is governed by the equation 

dt-Hri^) = {Ht,'Ht{^)} « I d^y{{vcC{y),nT{^)} + rc{C''{y),nT{^)}) = 
= [ dV{^cC(y),^T(x)} = I^C^^(x) « 

(3.41) 

using ( 3.37| ) and also using the fact that = 0. In other words 1-Lt is also preserved 
during the time evolution of the system without any restriction on the lapse function A^. 
Finally the requirement of the preservation of the constraint C^^ has the form 

dtC'\i^) = {Ht,C'\i^)] = 

d'y {{NnTiy),C''{^)} + vc{y) {C(y),C^^(x)}) = 

(3.42) 



using the fact that Fc = 0. Then with the help of the equation ( |3.40D we can argue that 
this solution can be solved for vc at least in principle. 

Let us outline our results. We have following first class constraints ttat ,7rj 
, "Hj ?a , y.T ~ together with the second class constraints C , C^^ ^ 0. Then we 
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also have 10 metric components gij,N, and corresponding conjugate momenta tt*-' , ttat , tTj 
and 4 scalar fields <f>^ with conjugate momenta pA- In general we have D = 28 phase space 
degrees of freedom. On the other hand we have S = 2 second class constraints and F = 8 
first class constraints. As a result the number of physical degrees of freedom is equal to 

m 

Np.d = {D-S-2F) = 10 (3.43) 

which is the correct number of the physical degrees of freedom of the massive gravity. In 
other words C and C^^ eliminate the ghost field and corresponding conjugate momenta at 
least in principle. 
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